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Abstract. We study four dimensional N = 2 supersymmetric gauge theory in the Q- 
background with the two dimensional 3\f = 2 super-Poincare invariance. We explain 
how this gauge theory provides the quantization of the classical integrable system 
underlying the moduli space of vacua of the ordinary four dimensional N = 2 theory. 
The e-parameter of the Sl-background is identified with the Planck constant, the 
twisted chiral ring maps to quantum Hamiltonians, the supersymmetric vacua are 
identified with Bethe states of quantum integrable systems. This four dimensional 
gauge theory in its low energy description has two dimensional twisted superpotential 
which becomes the Yang- Yang function of the integrable system. We present the 
thermodynamic-Bethe-ansatz like formulae for these functions and for the spectra 
of commuting Hamiltonians following the direct computation in gauge theory. The 
general construction is illustrated at the examples of the many-body systems, such 
as the periodic Toda chain, the elliptic Calogero-Moser system, and their relativistic 
versions, for which we present a complete characterization of the L 2 -spectrum. We 
very briefly discuss the quantization of Hitchin system. 
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1. INTRODUCTION 

It has been realized in the recent years [H [H O [H |5] that there exists an intimate 
connection between the vacua of the supersymmetric gauge theories and the quantum 
integrable systems. This connection is quite general and applies to gauge theories in 
various spacetime dimensions. In the short review articles [U[5] a large class of the two 
dimensional gauge theories were shown to correspond to the finite dimensional spin 
chains (and their various limits usually studied in literature on integrable models). 

We report here on the new developments - we establish the connection between the 
four dimensional supersymmetric gauge theories and the quantum many-body systems. 
More precisely, we consider the N = 2 supersymmetric gauge theories in four dimen- 
sions. We subject it to the O-background [6] in two out of four dimensions. The 
^-background is a particular background of the N = 2 supergravity in four dimen- 
sions. The general ^-background is characterized by two complex parameters £i,£2, 
which have the dimension of mass. These parameters were introduced in [U [7J and 
used in [8] to regularize the integrals over the instanton moduli spaces which arise in 
the supersymmetric gauge theories and the bound state problems in the supersym- 
metric quantum mechanics. It was suggested in [7] that by deforming the Donaldson 
supercharge Q to its equivariant version Q + VfGu the instanton partition functions 
would become computable (in fact, an example of an instanton integral was proposed 
in [7J) and could ultimately test the Seiberg-Witten solution [9] of the four dimensional 
IN" = 2 theory. This program was completed in [6] . It turns out that both the ideology 
and the specific examples of the integral formulae of [7J are important in our current 
developments related to Bethe ansatz of the quantum many-body systems. 

When both parameters E\ , £2 of the ^-background are non-zero the super-Poincare 
invariance is broken down to a superalgebra with two fermionic and two bosonic gen- 
erators corresponding to the rotations in R 4 . If one of these parameters vanishes, e.g. 
£2 = 0, then the resulting theory has a two dimensional 3sf = 2 super-Poincare invari- 
ance. This is the theory we study in the present paper. It is characterized by a single 
complex parameter, which we shall denote simply by e. 

Our main claim is: the supersymmetric vacua of this gauge theory are the eigenstates 
of the quantum integrable system obtained by the quantization of the classical algebraic 
integrable system underlying the geometry of the moduli space M v of undeformed N = 2 
theory. The Planck constant, the parameter of the quantization, is identified with the 
deformation parameter £. 

Recall that the four dimensional gauge theory with £\ = £2 = is characterized, 
at low energy, by the prepotential 3{a) which is known to be related to some classical 
algebraic integrable system. For example, the pure SU(N) N = 2 super- Yang-Mills cor- 
responds to the periodic Ajy_i Toda chain [10]. As shown in [11] the vector multiplet 
part M v of the moduli space of vacua, the Coulomb branch, of any N = 2 super- 
symmetric gauge theory, is a base of the Liouville fibration of some classical algebraic 
integrable system. Once the classical integrable system is identified one can write down 
the gauge theory low energy effective action. When £\ = £ ^ 0, £2 = the role of ^(a) 
in the gauge theory is played, in a certain sense, by the two dimensional twisted super- 
potential W(a;e). We show it also has a meaning in the algebraic integrable system, 



1 



NIKITA A. NEKRASOV AND SAMSON L. SHATASHVILI 



albeit in the quantum one. It is identified with the Yang- Yang counting function [12] 
governing the spectrum of the quantum system. One does not always know a priori how 
to quantize an algebraic integrable system, just knowing its classical version does not 
suffice. However, W(a; e) is computable by the gauge theory methods (and so is ^(a)), 
so the situation got reversed - the gauge theory helps to learn about the quantization 
of an integrable system. 

This is an important part of the general program, whose details will appear in the 
longer version [13] which in addition (and in particular) will combine the results re- 
ported here and those in [H [5] to form a unified picture. 

There are numerous applications of the correspondence [H El [31 HI [5]: the gauge 
theory applications, the study of quantum cohomology, (infinite-dimensional) repre- 
sentation theory, harmonic analysis, the many-body quantum mechanics. 

In a sense the most general non-relativistic algebraic integrable system is the so-called 
Hitchin system. Its quantization is an interesting and important problem, whose special 
cases, corresponding to the degenerate Riemann surfaces, are by now old and classical 
problems. In this paper we illustrate the power of our methods at the examples of 
two such degenerations, the quantum elliptic Calogero-Moser system, and its limit, the 
periodic Toda chain. We shall present a complete characterization of the L 2 -spectrum 
of these systems. 

The paper is organized as follows. The section 2 reviews the correspondence between 
the supersymmetric vacua of gauge theories with two dimensional N = 2 supersymme- 
try and Bethe vectors of some quantum integrable systems, and lays down a route to 
the analysis of the four dimensional theories. The section 3 is devoted to the four 
dimensional theories, the Q-background, the .F-terms, and the calculation of the effec- 
tive twisted superpotential. The section 4 reviews the integrability side of the story, 
first the classical algebraic integrable systems, then their quantization. The section 5 
describes the construction at the examples of the periodic Toda chain, elliptic Calogero- 
Moser system, the Ruijsenaars- Schneider model, and the Hitchin system. The section 
6 presents the thermodynamic Bethe ansatz-like formulae for W(a,e). 

Acknowledgments: The research of NN was supported by I'Agence Nationale de 
la Recherche under the grants ANR-06-BLAN-3.137168 and ANR-05-BLAN-0029-01, 
by the Russian Foundation for Basic Research through the grants RFFI 06-02-17382 
and NSh-8065.2006.2, that of SSh was supported by SFI grants 05/RFP/MAT0036 and 
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2. Review of the Bethe/gauge correspondence 

In this section we remind the correspondence between the supersymmetric vacua of 
the two dimensional N = 2 theories and the stationary (eigen) states of quantum inte- 
grable systems and then propose a four dimensional generalization. The correspondence 
can also be applied to the study of the topological field theory obtained by twisting the 
physical supersymmetric gauge theory. 

2.1. Twisted chiral ring and quantum integrability. The space of supersymmet- 
ric vacua of a theory with four supersymmetries carries a representation of a commu- 
tative associative algebra, the so-called (twisted) chiral ring, see e.g book [II] . 

For example, in two space-time dimensions, the 3Sf = 2 supersymmetry is generated 
by the fermionic charges Q±,Q±, which obey the anticommutation relations: 

{Q±,Q±} = 2(H±P), 

(2.1) {Q + ,Q_} = {Q + ,Q_} = 

{Q+,Q_} = {Q+,Q_} = 

the last two lines being valid in the absence of the central extension, induced, e.g. by 
some global symmetry charges. 

The twisted chiral ring is generated by the operators k = 1,2,... which (anti)- 
commute with the operator 

(2.2) Q A = Q+ + Q_, {Q A ,Q A } = H 

Analogously one defines the chiral ring, whose generators (anti)commute with the op- 
erator: 

(2.3) Qb = Q+ + Q-, {q b ,qU = h 

In our work we concentrate on the Q^-cohomology and assume that the possible central 
extension of (|2.ip leaves Q A nilpotent = 0. The local operators Ofc(x) are indepen- 
dent up to the QA-commutators of their location x. Their operator product expansion 
defines a commutative associative ring, 

(2.4) Q i Q j = <$ j G k + {Q A ,...} 

If |0) is a vacuum state of the Hamiltonian H, H\0) = 0, then so is Oj|0) = \i), and 
moreover the space of vacua is the representation of the twisted chiral ring. 

Thus the space of supersymmetric vacua, which can be effectively studied using 
the cohomology of the operator Q A (or Qb), is the space of states of some quantum 
integrable system: 

(2.5) Jjquantum = ker Q A / im Q A 

The actual vacua are the harmonic representatives of this cohomology as follows from 
(12. 2|) . (12, 3h . The operators and more generally the functions of O^'s are the quantum 
Hamiltonians. 
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The operators 0^ can be chosen under the assumption of the absence of massless 
charged matter helds to be the gauge invariant polynomials of the complex scalar a in 
the vector supermultiplet, 

One is looking for the common eigenstates of these Hamiltonians: 
(2.7) * A € H quantum , fc * A = E k (\)* x 

where -E'fc(A) are the corresponding eigenvalues, and A are some labels. In general they 
are complex, Ek(X) £ C. 

The important, or at least the interesting, problem is to identify the quantum inte- 
grable system given an N = 2 gauge theory, or to solve the converse problem - to find 
the N = 2 theory given a quantum integrable system. For a large class of models on 
both sides this problem has been solved in [U El El IH [5] . 

In most interesting cases the super symmetric gauge theory at low energies has an 
effective two dimensional abelian gauge theory description (with four supercharges), the 
so-called theory on the Coulomb branch. The supersymmetric vacua are determined 
in terms of the exactly calculable effective twisted superpotential W eS (a). Loosely 
speaking, given a vector of electric fluxes (m, . . . , n r ), with rij G Z the vacua are given 
by the critical points of the shifted superpotential W eS (a) — 2iri YJi=i n i° % (sometimes 
we loosely refer to the critical points of the superpotential without explicitly saying 
"shifted"): 

, . 1 dW cS (a) 

(2.8 — — = m 

v ; 2m da 1 

as follows from the consideration of the effective potential 

1 •• / dW° S \ ( dW eS 

(2.9) U H (a) = -g^ I -27rm 4 + ^— j +2vrin j + ^— 

Equivalently: 

, 2 ,o) ^(—L2j =1 

2.2. Topological field theory. The two dimensional N = 2 gauge theory can be 
topologically twisted. In the twisted version the supercharge Qa plays the role of 
the BRST operator. The vacua of the physical theory are the physical states of the 
topological theory. 

The action of the topological field theory can be brought to the simple form 

by the so-called quartet mechanism (one adds the anti- twisted superpotential iTra 2 , 
and sends i — > oo). Here F\ is a curvature of abelian gauge field A 1 and ip l is anti- 
commuting 1-form on S, the super-partner of A 1 . In this form the gauge theory becomes 
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a contour integral over the a field (the field a being eliminated). The original supersym- 
metric theory contains information hidden in the so-called Z)-terms, which ultimately 
leads to the wall-crossing phenomena [15j . 

The canonical quantization of the theory (|2.1ip on the cylinder S = R x S 1 is 
simple. Indeed, the only physical degree of freedom of this theory is the monodromy 
expzif = exp <|> sl A 1 of the gauge field around the circle S 1 and the momentum conju- 
gate to if, h = dW eS /d(Ji. Since if takes values in a circle (due to the large gauge 
transformations) the conjugate variable Ii is quantized, leading to the equations (12. 8|) . 



Our conclusion is that one can study the vacuum sector of the N = 2 gauge theory 
or one can study the topologically twisted version of the same theory - ultimately one 
deals with the Bethe states of some quantum integrable system. This correspondence 
benefits all three subjects involved. 

2.3. Yang- Yang function and quantum spectrum. In [21 [3] the form (|2.1ip of the 

effective action and the quantization argument above for the two dimensional theory 
studied in pQ was used to identify W cS (o-) with the Yang- Yang function of the N- 
particle Yang integrable system. However both the quantization argument and the 
identification of W eS with a Yang- Yang function are not restricted to this case only. 

It is a remarkable feature of the quantum integrable systems that their spectrum 
can be sometimes studied using Bethe ansatz. The spectrum of the quantum system is 
determined by the equations on the (quasi) momentum variables (rapidities) Aj, which 
enter the parametrization of the eigenfunctions. It is even more remarkable that the 
quasimomenta are determined by the equations which have a potential: 



The function Y(X) is called the counting, or YY function [T2], and (|2.12|) is the corre- 
sponding Bethe equation [IB] . 

It has been demonstrated in [H [S] that for many interesting gauge theories in two, 
three (compactified on S 1 ) or four (compactified on T 2 ) dimensions the equation (|2.10j) 
coincides with Bethe equation, for a large class of interesting quantum integrable sys- 
tems. Moreover, the effective twisted superpotential of the gauge theory equals the YY 
function of the quantum integrable system (e. g. for G = U{N)): 



when the Coulomb branch moduli Oi are identified with the spectral parameters Aj 
of the quantum integrable system. In addition, the expectation values of the twisted 
chiral ring operators in the vacuum |A) given by the solution a = A of (|2.8p coincide 
with the eigenvalues of quantum Hamiltonians of integrable system: 




1 dY(\) 
2tti d\i 



m e z 



(2.13) 




(2.14) 



(A|O fc |A) = £ fc (A) 



(2.15) H fc * A = £ fc (A)* A 

Most of these models of [H [5] are intrinsically quantum. Their Hilbert spaces are finite 
dimensional. The Planck constant in these cases is not a continuous variable. These 
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theories are not very natural from the point of view of the quantization procedure, 
where one is given a classical integrable system and is asked to construct the quantum 
integrable system. 

2.4. Quantization from four dimensions. In the current paper we study a novel 
type of theories, which originate in four dimensions. They have a continuous parameter 
e, which becomes the continuous Planck constant of a quantum integrable system. The 
corresponding quantum integrable systems have infinite-dimensional Hilbert spaces. 
For example, we shall give a solution to the quantum periodic Toda chain (pToda), an 
elliptic Calogero-Moser system (eCM) and their relativistic analogues. We shall also 
make some remarks on the general Hitchin system. 
Our strategy is the following: 

(1) Start with a four dimensional N = 2 gauge theory (for example, one may 
take a pure K = 2 theory with some gauge group G or the N = 2* theory, 
the theory with one massive hypermultiplet in the adjoint representation); its 
low energy effective Lagrangian is determined in terms of a single multi- valued 
analytic function of the Coulomb moduli (a\, . . . ,a r ), r = rank(G), called the 
prepotential, 5"(a; m, r). Here m, r etc. are the parameteres of the gauge theory, 
whose meaning will be clarified later on. 

(2) The theory in the ultra-violet has the observables 0^ 

where <f> is the complex scalar in the vector multiplet. The observables 0^ cor- 
respond to the holomorphic functions on M v . These observables are singled 
out by their (anti-)commutation with the supercharge Qa (this supercharge be- 
comes the Q-operator of the Donaldson- Witten theory in the standard N = 2 
twist). One can formally deform the theory in the ultra-violet, 

(2.16) J tree -» J tree + tkQk 

k 

In the low energy the deformed theory is given by the family of prepotentials 
3(a\t;m, r), t = (t*.). For some models the deformed prepotentials were com- 
puted in [151 El HHJ- For our purposes the deformation can be studied 
formally, i.e. without addressing the convergence issues. However, some of the 
couplings in t can be interpreted as shifting the ultraviolet complexified gauge 
coupling r. For these couplings a finite deformation can be studied. For other 
couplings only the first order deformation is needed (and therefore the con- 
tact term problem [15] does not arise). The deformed prepotential 5F(a;t;m, r) 
generates the vevs of O^'s: 

m k " {0k)a 

(3) The prepotential 9~(a, m, r) has an interpretation in terms of a classical alge- 
braic integrable system where it plays the role of the generating function of a 
Lagrangian submanifold L of a symplectic complex vector space relating two 
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types of complex action variables, a 1 and a^i- The t-couplings deform L to 
another Lagrangian submanifold Lt [15j . 
(4) Finally, we introduce one more deformation parameter, e, which corresponds 
to subjecting the theory to the ^-background S 1 x R 1 x R £ 2 , or R 2 x (see 
the next subsection). The F-terms of the low energy Lagrangian are now two 
dimensional with the twisted superpotential which we denote by W(a; t; e). For 
the discussion of the supersymmetric vacua only these twisted F-terms are 
relevant, therefore the effective description of the low energy physics is two 
dimensional. In particular, one derives the equation determining the vacua, as 

(2.17) V = 2vrini , m € Z 

while the spectrum of the twisted chiral observables is computed from: 

0W(a;t;e) 



(2.18) 



dt k 

t=o 

We shall note that the equation (|2.17p can be written, in the examples studied in 
this paper and probably in more general situations, in terms of the factorized 
5-matrix of the associated "hyperbolic" many body system corrected by the 
"finite-size" terms (i.e. the instanton corrections in q in the case of N = 2* 
theory, the corrections in A for the pure Ji = 2 theory etc.): 

(2.19) — - — - — = i h y log S(a,k — a 7 ) + the finite size terms 

oa k e 

where we set t k = for k > 2. 
(5) As in [H El El H [5] the equations (l2Tfl) . (l2~18l) determine the spectrum of 
the quantum integrable system, where e plays the role of the (complexified) 
Planck constant. We identify this system with the quantization of the classical 
algebraic integrable system describing the low energy effective theory of the 
four dimensional 3Sf = 2 theory. 
Remark. Note the unfortunate notational conflict. The order parameters, <Tj, the 
eigenvalues of the complex scalar in the vector multiplet, are denoted traditionally 
by G{ in the context of two dimensional gauge theories, by a 1 in the context of four 
dimensional gauge theories, by (pi in the context of topological gauge theories. The 
same parameters exhibit themselves as the Bethe roots in our correspondence with the 
quantum integrable systems. In that world they are denoted by Aj. In the context of 
the periodic Toda chain these variables are denoted by <5j in [19] , and by t a in [20J . We 
failed to propose a unified notation and followed the traditional way of denoting the 
same object by different letters corresponding to the different context. We hope the 
reader will not be too confused by this. 
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3. Four dimensional gauge theory 

In this section we describe in some detail the fi-deformation of the supersymmetric 
gauge theory. We analyze the effective theory and observe that it can be related, at 
the level of cohomology of some supercharge, to a supersymmetric gauge theory in two 
dimensions. The two dimensional theory has an effective twisted superpotential which 
we analyze. 

3.1. The O-background and twisted masses. A quantum field theory in k + 2 
space-time dimensions can be viewed, formally, as a two dimensional theory with an 
infinite number of fields. If this theory is studied on flat k + 2 dimensional space-time, 
or on a space-time fibered over a two-dimensional manifold £ with the Euclidean fibers 
R fc , it has (in addition to the global symmetries of the k + 2-dimensional theory) a 
global symmetry group E(&) of isometries of the /c-dimensional Euclidean space R fc . 
Accordingly, if the theory has anN=2 supersymmetry on S, then one may deform it 
by turning on the twisted masses corresponding to the global symmetry E(k). Unlike 
the conventional global symmetries which typically form a compact Lie group with 
the unique, up to a conjugation, maximal torus, the group E(&) has several, for k > 
1, inequivalent Cartan subgroups. Thus there exists several physically inequivalent 
deformations of the k + 2 dimensional theory. 

For example, one may choose a subgroup R fc_2 ' x SO{2) 1 of translations in H k ~ 21 and 
rotations in the I orthogonal two-planes. The theory with twisted masses corresponding 
to the translations in R fc_2 ' is equivalent to the ordinary Kaluza-Klein compactification 
on a torus T fc ~ 2 '. Such theories are studied in [H [5]. It is shown there, that if one 
starts with the four dimensional gauge theory with N = 2 supersymmetry, with L 
hypermultiplets in the fundamental representation, and the gauge group U(N), then, 
upon the compactification on the two-torus T 2 one gets, via (|2.5p the su(2) XYZ spin 
chain, where the elliptic parameter of the spin chain is identified with the complex 
structure modulus r of the compactification torus. The three dimensional theory gives 
rise to the XXZ spin chain, and the two dimensional one corresponds to the XXX spin 
chain. The number ./V of colors is equal to the excitation level of the spin chain (the 
total spin's projection S z is equal to N — while the number of flavors L is the length 
of the spin chain. 

3.1.1. Four dimensional theory on R 2 . Another possibility is to consider the twisted 
masses corresponding to the rotational symmetry. In this case one gets the theory in 
the ^-background [6]. 

Consider a four dimensional 3sf = 2 theory on a four manifold M 4 fibered over 
a two dimensional base £ with the f2-background along the fibers R 2 . Somewhat 
schematically we shall denote the fibers by R 2 . The base X of the fibration could be a 
two-plane R 1,1 or a cylinder S 1 x R 1 . One can also study the twisted theory for which 
£ could be an arbitrary Riemann surface. 

In our main examples, corresponding to the periodic Toda and elliptic Calogero- 
Moser systems our staring point would be the pure !N = 2 super- Yang-Mills theory, or 
the 2NT = 2* theory, corresponding to the gauge theory with a single adjoint massive 
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hypermultiplet. In the limit of vanishing mass of the adjoint matter fields the latter 
becomes an N = 4 theory. 

Now let us give some details on the fi-deformation. Let us denote the coordinates 
on the fiber R 2 by (x 2 ,x 3 ), and the coordinates on the base S by (x ^ 1 ). Introduce 
the vector field 

(3.1) U = x 2 d 3 - x 3 d 2 

generating the U(l) rotation in R 2 . Let e G C be a complex parameter and let 
V = eU, V = eU be the complex vector fields on R 2 . 

The bosonic part of the pure 3sf = 2 super- Yang-Mills Lagrangian of the theory on 
R 2 is simply: 

(3.2) L = — 3_tr F A *F + tr {D A (j) - et v F) A * (D A <j) - eluF) + 

+itr ([0, 0} + l v D a [e4> -ecf}) 2 + ^-txF A F 

It is clear that the Poincare invariance in the (x ^ 1 ) directions is unbroken, and it is 
possible to show that in fact the two dimensional N = 2 super-Poincare invariance is 
preserved. Thus there are four supercharges. 

The only (twisted) F-terms of the low energy effective theory are two-dimensional 
and can be represented as the non-trivial twisted superpotential W(a; e) (where as 
before a denotes the complex scalar in abelian vector multiplet). 

If we send e back to zero then the low energy theory is fully four-dimensional with 
a continuous moduli space M v of vacua and the low energy effective Lagrangian is 
described in terms of the prepotential 7 {a). 

For non-zero e the theory has a discrete set of vacua given by the minima of the 
potential (|2.9p . (|2.10p . These vacua are therefore the solutions to the equation ([2.170 . 
which we shall later on identify with the Bethe equation of some quantum integrable 
system. 

3.1.2. Calculation of the twisted superpotential. Our discussion would have had a rather 
limited significance were it not for the possibility of exact computation of W(a;e). 

Let us briefly explain our strategy. Consider the four dimensional theory in the gen- 
eral fi-background, with both rotation parameters £\, £2 non-zero. Then the effective 
theory has a prepotential 3~(a, £1,62) which is analytic in £1,62 near zero and becomes 
exactly the prepotential of the low energy effective four dimenional theory in the limit 
£i,e 2 -> 0: 

(3.3) Sf d = J ?W(a,e 1 ,e 2 ) + {Q,...} 

where we denote collectively by a all the vector multiplet scalars as well as background 
scalars, such as the masses of matter fields. At the same time, had we started with a 
two dimensional theory which is characterized by some twisted superpotential W(a), 
and had we subjected it to the two dimensional ^-background with the parameter £2, 
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the effective action would have had the form: 

(3.4) Sf d = J W^(a,e 2 ) + {Q,...} 

where W(a,£ 2 ) ~~ * W(a) as £ 2 — > 0, and Q is a certain supercharge which we use to 
study the vacuum states of our theory. In the Eqs. (13. 3h . (|3.4p the notations i* 1 ^, 
Vl-^ 2 ) refer to the cohomological descendents of the local operators F, W, etc. Now the 
standard manipulations with equivariant cohomology give: 

(3.5) f ^ 4 \a,e u e 2 ) = -[ & 2 \a,e 1 ,e 2 ) = J_y( a ,e 1} £ 2 ) 

7R4 £l JR2 £l£2 

modulo Q-exact terms at each step. By carefully manipulating the Q-exact terms we 
can connect the computation of the gauge theory partition function in the ultraviolet, 
which is given by a one-loop perturbative and a series of exact instanton corrections, 
to the computation in the infrared, using the Wilsonian effective action, where the 
effective energy scale can be sent all the way to zero: 

(3.6) exp— !— g r (a,£i,£ 2 ;q) = Z(a,£i,£ 2 ;q) = % pert (a, £i, £ 2 ; q) x 2, inst (a, £ X , £ 2 ; q) 

£l£2 

where we restored the dependence on the complexified bare gauge coupling 

•& 4iri 

q = exp 2-kit t = - — | 5- 

2vr g 2 

Now, by comparing the Eqs. (j3.3|) . (|3.4j) . (j3.6h we conclude: 

(3.7) W(a;£;q) =Limit £2 ^ [e 2 logZ(a, e 1 = e, £ 2 ; q)] = W pcrt (a; e; q) + W inst (a; e; q) 
Here the instanton part has an expansion in the powers of q, 

oo 

W inst (a;£;q) = q fc Wi nst (a; £) 
k=l 

while the perturbative part has a tree level term, proportional to log(q) and the one- 
loop term, which is q-independent. 

It follows that in the limit £ — > the twisted superpotential W(a, e) behaves as: 

(3.8) W(a;£;q) = ^^ + ... 

with . . . denoting the regular in e terms, and now the equations on the supersymmetric 
vacua assumes the Bethe (I2.17P form with the superpotential W(a,e) (13. 7D . (13. 8p . 

Our course is now pretty much set. We shall use the gauge theory knowledge of the 
instanton partition functions Z(a, £i,£ 2 ;q) to extract, via (13. 7j) . the twisted superpo- 
tential W(a, e;q). The details of this procedure are reviewed in the Section 6 where 
various ways of writing Z(a, £i,£ 2 ;q) and extracting W(a,£;q) are presented. For the 
purposes of the current discussion, the function W(a, e; q) is known and the logic above 
gives the desired equation for the supersymmetric vacua in the form of Bethe equations 

By the philosophy of [H EJ EJ [U [5] this superpotential serves as the Yang- Yang 
function of some quantum integrable system. In the present context the parameter e 
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plays the role of the Planck constant, and can be tuned to zero. In this limit we shall 
be able to use the quasiclassical asymptotics (I3.8P to identify the classical integrable 
system, whose quantization (in e) is the quantum integrable system in question. 

To this end we need to remind the role of the prepotential &{a; q) in the world of 
classical integrability |21j . 

Remark 1. Note that in the derivation [22] of the Seiberg-Witten prepotential 3~(a) 
from the direct instanton counting one evaluates the small £i,£2 — * asymptotics of 
Z(a, e±, £2; q) by a discrete version of the saddle point method, which connects nicely 
the theory of instanton integrals to the theory of limit shapes and random geometries. 
In our story we need to go beyond that analysis, see Section 6. Our results suggest 
that the "quantum theory of the limit shape" is related to the thermodynamic Bethe 
ansatz [HI [231 [Ml EHl EE [27] . 

Remark 2. The simplest case of the ^-background is in two dimensions. The two 
dimensional fi-background is characterized by the single parameter e. The partition 
function I a p{e) depends on two (discrete) parameters: the choice of the boundary 
condition at infinity, and the choice of the twisted chiral ring operator a inserted at 
the origin. These partition functions were studied in the context of the (equivariant) 
Gromov-Witten theory [28], where the ^-dependence comes from the coupling to the 
two dimensional topological gravity. When the theory is deformed in a way analogous 
to (|2.16p . the partition function becomes a matrix- valued function I a p{t; e) which solves 
the quantum differential equation: 

(3-9) ^(t; £ ) = C;(t)^(t; £ ) 

where C\ are the structure constants of the twisted chiral ring. 
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4. Integrable SYSTEMS 

In this section we remind a few relevant notions in the theory of classical and quantum 
integrable systems and in particular explain the role of prepotential in classical algebraic 
integrable system. We introduce the main examples, the periodic Toda, the elliptic 
Calogero-Moser systems, their relativistic versions, and the Hitchin system. We also 
introduce (formally) the t-deformations of these systems, along the lines of |15j . 

4.1. The classical story. The classical Hamiltonian integrable system is the collec- 
tion C?,lo, H), where IP is a 2n-dimensional smooth manifold endowed with the non- 
degenerate closed two-form u> and a collection H = (H\,H2, ■ ■ ■ ,H n ) of (generically) 
functionally independent functions H : IP — ► R n , which mutually Poisson-commute: 
{Hi, Hj} = 0, % ^ j. Here {A, B} = (w -1 ) d^Ad u B. We can view IP as a Lagrangian 
fibration H^^UcR" 

The classical Liouville- Arnold theorem states that if the common level set H -1 (/i) is 
compact, then it is diffeomorphic to the n-dimensional torus T n . Moreover, if H —1 (7i) 
is compact for any h in a neighborhood U of a point ho € R n , then H~ 1 (U) is sym- 
plectomorphic to the neighborhood of a zero section in T*T n . One can then find the 
special Darboux coordinates (I, 92), I = (Ji, 1%, . . . , I n ), tp = (ip±, . . . , ip n ), called the 
action-angle variables, s.t. the Hamiltonians Hi, i = 1, ... ,n, depend only on I, Hi(I), 
while ipi are the periodic angular coordinates on T n with the period 2ir. Explicitly, the 
action variables are given by the periods 

(4.1) h = hi pdq 

J Ai 

of the one-form pdq = d -1 u; (one can give a more invariant definition), over some 
Z-basis in #i(T n ,Z). 

The notion of the classical real integrable system has an interesting complex analogue, 
sometimes known as the algebraic integrable system. The data (IP, to, H) now consists 
of the complex manifold IP, the holomorphic non-degenerate closed (2, 0) form lo, and 
the holomorphic map H : IP — ► C n whose fibers 3h = H _1 (/i) are Lagrangian polarized 
abelian varieties. The polarization is a Kahler form w, whose restriction on each fiber 
is an integral class [w] G H 2 (Jh, Z)nH 1,1 (Jf l ). The image IB = H(IP) is an open domain 
in C n . It has a special Kahler geometry, with the metric 

1 n 

(4.2) ds 2 = - Vim (da l ®da D A 

7T z — ' 

i=l 

where the special coordinates a l ,ao,i are given by the periods: 

(4.3) a % = — ([ pdq, am = — f pdq 



27r J A, " 27r JB, 



over the A and -B-cycles, which are the Lagrangian (with respect to the intersection 
form given by [w]) subspaces in Hi(Jh,Z). It follows that the two-form Y2i A 
da_D,i vanishes on H> thereby embedding the covering II of the complement U3\£ to the 
discriminant S C IB of the singular fibers to the first cohomlogy i? 1 (J/ l0 , C) of the fiber 
over some distinguished point /iq € IB, as a Lagrangian submanifold L. As such, it 
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comes with the function 3~ : £ —* C which can be locally viewed as a function of a 1 , 
such that 

(4.4) a Cjl = ^ 

The comparison of the Eqs. (|4.3p and (|4.ip suggests both a 1 and a^j are the complex 
action variables. Since the 2n-dimensional symplectic manifold has at most n function- 
ally independent Poisson-commuting functions, there ought to be a relation between 
a 1 and a^^'s. It is remarkable that this relation has a potential function. The action 
variables a 1 come with the corresponding angle variables = en + TijP J , Tjj = 9^3", 

while a£),j correspond to 4 >D,% = (t -1 )*"' 

(4.5) u = ^da i Ad^ = ^da Ai Ad^- ! 

i i 

Here aj,/3 l G R/27rZ are the real angular coordinates on the Liouville torus. 

4.1.1. The t- deformation. An algebraic classical integrable system (T, w,H) can be 
deformed in the following way. Consider a family (!Pt,Wt 5 Ht) of complex symplectic 
manifolds with the Lagrangian fibration given by the "Hamiltonians" Ht : CPt — * C n , 
over a (formal) multidimensional disk T>, parameterized by t. Then the variation of 
the symplectic form ut in t can be described by the equation: 

(4.6) J-„ t = Ly^ 
where 

.R 



doti A d/3 l + ^2 Imnjda 1 A da j 



i=l ij 

is the Kahler (1, l)-form, and 14 is the holomorphic Hamiltonian (in the uit symplectic 
structure) vector field corresponding to the Hamiltonian of the original integrable 
system. 

4.1.2. Quantization. The quantization of the classical integrable system is a (possibly 
discrete) family (yi e ,!H e ,H), of the associative algebras A £ , which deform the algebra 
of functions on the Poisson manifold (X, a; -1 ), the (Hilbert) vector spaces "K e , with the 
action of A e , and the operators H = (Hi, . . . , H n ), Hi GA £ , which mutually commute 
[Hi,Hj] = 0, i / j, and generate "K e in the following sense: the common spectral 
problem 

(4.7) Htf> = 

defines a basis in "K e . Here e has the meaning of the Planck constant. 

The construction of the common eigenstates and the spectrum of the operators &% is a 
problem of the coordinate Bethe ansatz, quantum inverse scattering method (Algebraic 
Bethe ansatz [3Q |. 131 } 132]). quantum separation of variables [33], Baxter equation [34], 
the approach of [35J to the spectral curve quantization and various other versions of 
the Bethe ansatz. 
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In this paper we take a different route - via the supersymmetric gauge theory, along 
the lines of [B El SJ [5] . The gauge theory allows to find the exact spectrum of (|4.7j) 
which is the invariant of the choice of polarization used in the quantization procedure. 
The closest to our approach in the integrability literature seems to be that of [35j . 

In our story the algebra A e is the deformation of the algebra of holomorphic functions 
on y. We shall assume the existence of the global coordinates Pi,Xi (we do not assume 
them to be the globally defined holomorphic functions on CP, but on some covering 
space). The algebra A e is generated by pi,Xi, obeying 

(4.8) \Pi,Xj] = e5ij , 

"K £ is defined as the space of appropriate holomorphic functions of x«, and the repre- 
sentation of A £ is given by: 

d 

(4.9) Xi = Xi, pi = e— 

OXi 

The next question is the construction of the Hilbert space where A e is represented. 
If y were a cotangent bundle to a complex manifold M, the algebra A £ would be 
isomorphic to the algebra of holomorphic differential operators on M. However, it is 
rarely the case that there are interesting holomorphic differential operators which are 
defined globally on M. Moreover, the naive complexification of the quantization of 
the quantization of T*Mr (which produces the differential operators on Mr acting 
in the space of half-densities L 2 (Mr; K^jf)) produces the K l J { 2 - twisted differential 

1/2 

operators. It may well happen that the space of global sections of the KJ. - twisted 

1/2 

differential operators is non-zero yet the space of global sections of KJj where these 
operators would have acted is empty. This is the situation with the quantization of 
Hitchin system as discussed in [36j . 

In our story the algebra A £ is the noncommutative deformation of the algebra of 
holomorphic functions on 7, yet it is represented in the regular L 2 -sections of some 
line bundle on a real middle dimensional submanifold IPr of 7. The choice of !Pr is 
apparently made by the boundary conditions in the gauge theory. We do not have a 
complete understanding of this issue yet, but let us make an important: 

Remark: The equation (|2.8p and the discusion in Section 2 show that in our 
approach we quantize, for the type A model (defined below) the real submanifold 
which projects onto the locus Re (<9W(a; t; e)/da l ) = in the base, and cuts out a 
middle dimensional real torus in the Liouville fiber; for the type B model it projects 
onto the locus Im (a 1 /e) = 0. 



4.1.3. Quantization and t- deformation. The quantum integrable system can be de- 
formed by making Hi depend on the additional parameters t = (ti, . . . , t n ), in parallel 
with the classical deformation (|4.6p . so that they define a flat connection depending on 
a spectral parameter k: 

(4.10) [ K Hi (t), K - Hj {t)}=0 
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for all i,j = 1, . . . ,n. This deformation, when applied to quantum Hitchin system, is 
related to the KZB connection in the WZW conformal field theory. The gauge theoretic 
meaning of the ^-parameter is 82 of the general fi-background. 
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5. Examples 
We now proceed with explicit examples. 

5.1. The periodic Toda chain. 

5.1.1. The classical system. The periodic Toda chain is the system of N particles 
xi,. . . ,xn on the real line interacting with the potential: 

(N-l 
£ e^-^+i + e XN ~ Xl 
i=i 

The phase space of this model is IPr, = T*H N , with the coordinates (pi,Xi)f =1 , where 
Pi,Xi £ R, the symplectic form uj = J2iLi ^Pi A dxj, and the Hamiltonians 

(5.2) Hi = ^ Pi 

i 

(5.3) H 2 = l -Y J Pf + U{xu...,x N ) 
(5.4) 

(5-5) H k = + - 

(5.6) 

The complexified Toda chain has the phase space 7 = T* (C X ) N , with the coordinates 
(pi,Xi)f =1 where pi GC,x,G C/{2m)Z. 

To describe this model as the algebraic integrable system we introduce the Lax 
operator 



/ pi A 2 e Xl_Xa e~ z \ 

1 p 2 A 2 e X2 ~ X3 

1 ps AV 3- * 4 o 

(5.7) $(z) = ... ... 

... ... 

... ... ... pjv-i A 2 e XN - 1 ~ XN 



\K 2 e XN ~ xl e z ... 1 p N J 

and define the Hamiltonians hi, . . . , h^ as the coefficients of the characteristic polyno- 
mial: 

(5.8) Det (x - = -A 2 V - e~ z + x N + fuz*" 1 + fore*" 2 + . . . + h N 

The Hamiltonians /ii , hi are then given by 

N 

hi = Pi 
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and 

(5.9) h 2 = - y~]piPj + U(%i, . . . ,x N ) 

i<j 

which differs from H2 defined in f)5.3f) by the term \h\. 

Define the spectral curve C C x C x , i£C,z£ C/2iri, as the zero locus of the 
characteristic polynomial (15.8|) . For each value h = (Hi, H2, ■ ■ ■ , Hjsr) this is a curve 
which is a genus N — 1 hyperelliptic curve with two points where x = 00 deleted. The 
fiber H _1 (/i) is given by the product Cx^h- The C-factor corresponds to the center-of- 
mass mode J2% x «i the compact factor Sh = Jac(C/ l ) is the Jacobian of the compactified 
curve 6^. 

The complex action variables a 1 , ao,i can be computed as the periods of the differ- 
ential 

(5.10) A = —xdz 

5.1.2. The gauge theory. The gauge theory significance of the periodic Toda chain is 
that the potential ^(a) in the Eq. (|4.4p defined using the family of spectral curves (|5.8p 
and the differential (|5.10j) coincides with the prepotential of the low energy effective 
Lagrangian of the pure N = 2 gauge theory with the gauge group U(N). 

5.1.3. The quantum system. The quantization of the periodic Toda chain is achieved 
by promoting h^s defined by (|5.8p to the differential operators acting on functions 
of (x\, . . . ,xn), Pi = zd Xt . It is possible to show that the potential normal ordering 
ambiguities do not arise in this case. One is looking for the eigenfunctions of the form: 

(5.11) *(a?i, • • • , x N ) = e x ipk(xi - x,x 2 - x, . . . ,x N - x) 
where 

1 N 
i=i 

tpk € L 2 (R N ~~ 1 ). When e = —ih and A,fi£R this problem has a discrete real spectrum 
(for fixed k). We are after the effective characterization of this spectrum. It turns out 
that this model allows a complex analytic continuation in the parameters e, A, so 
that the spectrum remains discrete, yet in general complex. We shall call this spectral 
problem the type A quantum periodic Toda. 

The quantum periodic Toda admits another, somewhat unconventional (for N > 2) 
formulation, which also leads to the discrete yet complex spectrum. In this formulation 
we make the differential operators act on functions ^(x\, ... ,xjsr) which are 2iri- 
periodic, and non-singular for some fixed value of Rexi, . . . ,HexN- We shall call this 
spectral problem the type B quantum periodic Toda. 

Note that for N = 2 case the type B periodic Toda is equivalent to finding the (quasi)- 
periodic solutions of the canonical Mathieu's differential equation, while the type A 
model corresponds to the L 2 solutions of Mathieu's modified differential equation. 

Remark. The quasi-periodic solutions of the differential equations are parametrized 
by the so-called Bloch-Floquet multipliers, which in our correspondence come from the 
two-dimensional theta angle which might come from the 5-field in four dimensions, 
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or from the peculiar deformation of the four dimensional tree level prepotential cx 
e$tr$, in the O-background H 2 . Of course we can set i? = and discuss the periodic 
wavefunctions. This remark applies to all the many-body systems. 

5.2. Elliptic Calogero-Moser system. 

5.2.1. The classical system. The elliptic Calogero-Moser system (eCM) is the system 
of N particles xi,x 2 , ■ ■ ■ ,xjy on the circle of circumference j3, i.e. x% ~ + /?, which 
interact with the pair-wise potential 

(5.12) U(xi,X2, ■ ■ ■ ,xn) = m 2 '^2u(xi - Xj) , 

u(x) = C(J3) + V ^ = -dl log e(x) 

V ' V 1 ' ^sinh 2 (x + fe/3) x 6 V ' 

where is the odd theta function on the elliptic curve E T with the modular parameter 

n ' 

e(x) = - V (-l) k q^e 2kx , q = exp27rir 



feez 



2 



and C(j3) is some constant which depends on 0. 

Again, there exists a Lax representation [37J of the eCM system: 

(5.i3) « B w=M» + m efalx^w {1 - S,l) 

In the limit f3 — > oo, m — > oo, such that A 2N = m 2N q is kept finite, the eCM system 
becomes the periodic Toda chain [38], where 

x fM = _^ + x pToda 

5.2.2. The gauge theory. The spectral curve Det(<J>(z) — x) = 0, which is an iV-sheeted 
ramified cover of the elliptic curve E T where z lives, is the Seiberg-Witten curve of a 
remarkable four dimensional N = 2 theory. This is an SU(N) gauge theory with a 
massive adjoint hypermultiplet. The parameter m in (15, 13ft is the mass of the adjoint 
hypermultiplet, the complex structure of the curve E T is determined by the complexified 
bare gauge coupling of the ultraviolet theory (which is in fact the superconformal "N = 4 
super- Yang-Miills). We shall call this theory the N = 2* theory. 



5.2.3. The quantum system. 

e 2 N d 2 

(5.14) H 2 = - m{m + e)^2 p(xi - x s ) 

i=l 1 i<j 

When s = —ih, m = —ve, and h,u,/3 £ R+ the spectral problem of the operator 
(|5.14p is well-known: one is looking for the (quasi-)periodic (in (5) symmetric functions 
*&(xi, . . . ,xjsr), which are non-singular in the fundamental domain, and vanish at the 
diagonals: 

(5.15) *(xi, . . . ,x N ) ~ (xi - Xj) v 
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It follows from our results that this problem has an analytic continuation where e, m, (3 
all become complex. The monodromy of this continuation is an extremely interesting 
problem which we shall not be able to cover in this exposition. 

Just like in the case of the periodic Toda chain one may consider various spectral 
problems for the Hamiltonian Hi and its higher order counterparts. The pleasant bonus 
of having an elliptic potential is the similarity of these problems. The type A quantum 
elliptic Calogero-Moser is, therefore, a problem which we just described: finding the 
/3-periodic I? (on the subspace with fixed Im(xj//3)) wavef unctions, with the (15. 15H 
behavior. The type B problem is a problem of finding the iri- periodic I? wavefunctions 
(on the subspace with fixed Re(xj)). The SL2CZ1) transformation r — > — £ maps the 
type A problem to the type B problem and vice versa. 

Superficially, however, one may worry that the type A and the type B problems are 
quite different in nature. In the limit (5 — > 00 the type A problem becomes that of the 
hyperbolic Sutherland model, which has a continuous spectrum. In the same limit the 
type B problem becomes that of the trigonometric Sutherland, which is well-studied, 
has a discrete spectrum, and Jack polynomials (up to a ground state factor) as the 
eigenf unctions. Our approach connects all these models. 

In the N = 2 case we are solving the celebrated Lame's differential equation |39j. 

5.3. Hitchin system. The previous systems are the degenerate examples of the so- 
called Hitchin integrable system. Its phase space 7 is the partial resolution of the 
cotangent bundle to the moduli space M of semistable holomorphic bundles on a com- 
plex curve C. More precisely, the phase space is the space of pairs (Vg,<& = & z dz), 
where = 8 + A, A = A 2 dz, is the (0, l)-part of a connection on a vector bundle E 
over C, and 5> 2 is the Higgs field, a (l,0)-form valued in the endomorphisms of the E. 
The vector bundle E gets a holomorphic bundle structure by declaring the solutions to 
VjX = to be the holomorphic sections. The Higgs field Q z must obey: 



Divide by the complex gauge transformations, or make a symplectic quotient with 
respect to the compact gauge transformations. In the latter case one imposes the real 
moment map condition: 



The Hamiltonians are obtained as follows: consider the (j, 0)-differentials on C, given 
by tr$ J . Due to (|5.16p these are holomorphic j-differentials, and (for j > 1) there are 
(2j — 1)(<7 — 1) linearly independent (over C) such differentials. The base 23 of the 
corresponding integrable system is the vector space: 



(5.16) 



= 0, i.e. d- z $ z + [Ag, *J = 



(5.17) F A + [$,$] = 

The holomorphic symplectic structure is induced from 



(5.18) 





while the fiber over a point h G H is the Jacobian of the spectral curve: 
(5.20) C C T*C : = Det ($ - x) 
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(where we view x as the canonical Liouville one- form on T*C). The homology class [6] 
spectral curve is equal to iV[C]. It follows that the self-intersection number 2g(Q) — 2 = 
C.e = N 2 C.C = N 2 (2g - 2), therefore the genus 

(5.21) g(e) = l + N 2 (g-l) 

is equal to the dimension of 25. The polarization comes from the Kahler form: 

(5.22) zu= f tr (6 A A 5 A + 5<S> A <5$) 

Jc 

Let us now explain that the systems we considered so far are the degenerate cases of 
Hitchin system. Imagine we study a Hitchin system for the group U(N) on a curve of 
genus two. Now let us degenerate the curve, in such a way that it becomes a union of 
two elliptic curves, connected by a long neck. Then the equations (|5.16p can be solved 
on both elliptic curves independently, the only memory of the original curve being a 
boundary condition at the puncture representing the neck. The limiting Hitchin system 
would actually split as a union of invariant submanifolds, labelled by the coadjoint 
orbits of the complexified gauge group, attached to the puncture (as a hyper kahler 
manifold it fibers over (t <g> R 3 )/W) gO]. The limiting equations (EU6l) look like: 

(5.23) d- z $ z + [A- z ,$ z ] = J5 { V{z,z) 

where J represents a (complex) moment map of the group action on a coadjoint orbit 
0. Now let us consider a very special vase, where the orbit is diffeomorphic to 
T*CP Ar_1 . The corresponding moment map J is the N x N complex traceless matrix 
of which N — 1 eigenvalues coincide, 

J = mdiag (N - 1, -1, . . . , -1) . 

By solving the equation (|5.23j) one gets (up to an irrelevant gauge transformation) the 
Lax operator (|5. 13j) in the gauge where 

A 2 = diag (x 1 , . . .,x N ) 

is a constant diagonal matrix |41j . 

5.3.1. The quantum system and the gauge theory. We do not have much to say about 
these two topics in the case of Hitchin system for compact Riemann surface. The 
supersymmetric field theory realization of this system involves the six dimensional (2, 0) 
theory compactified on the Riemann surface C. Using the results of [32] one may hope 
to formulate the result as the four dimensional "N = 2 supersymmetric gauge theory, 
at least for the curves C which are close to the degeneration locus. This gauge theory 
can be further subject to the f2-background. As a result one would get a quantization 
of Hitchin system, together with an expression for its Yang- Yang function. It would be 
interesting to understand the relation of our approach to the constructions of [33] . 

5.4. Relativistic systems. So far we considered the systems which have a Galilean 
symmetry group, the space translations are generated by H± while the time translations 
are generated by H2, and there is a generator of the the Galilean boost S, such that 
[S, H2] = H\. The point is that the Hamiltonians H\, H2, ■ ■ ■ , Hjy are polynomials in 
the momenta pi, ■ ■ ■ ,pn- 
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There exists |444 H5] a one-parametric deformation of the eCM integrable system, 
whose Hamiltonians are trigonometric (hyperbolic) functions of the momenta. In par- 
ticular, the Hamiltonians P = Hi,E = H2 become 

N 



The quantization of the system (|5.24|) produces a family of commuting difference opera- 
tors. These operators, in the trigonometric limit, reduce (after a similarity transforma- 
tion), to the famous Macdonald difference operators [3B], whose eigenfunctions are the 
so-called Macdonald polynomials. In the context of a physical problem, the polynomial 
eigenfunctions correspond to the compact model, where the particles xx,...,xn live 
on a circle. The problem of interest for us reduces, in the trigonometric limit, to the 
hyperbolic system where the particles live on a real line, and the spectrum, as in the 
non-relativistic case, becomes continuous. As far as we know, there is no satisfactory 
treatment of this model in the literature. Therefore our results, presented below, might 
be of additional interest for the community working in the domain of harmonic analysis. 

In our analysis, the relativistic systems (|5.24p correspond to the five dimensional 
gauge theory compactified on a circle [47], subject further to the O-background in R|. 



(5.24) 




i=i 



N 



E = ^cosh(/3pi)/i(x) , 
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6. Superpotential/ Yang- Yang function W(a,e;q) 

In this section we show how W(a, e; q) is computed and give several representations 
for it. The W(a, e; q) from this section needs to be inserted in (12,17h . (|2.18p in order to 
write the explicit expression for the spectrum. 

6.1. Thermodynamic Bethe Ansatz. The perturbative part of W is written using 
the special function w £ (x). It obeys 

—w e {x)=\ogT[l + - 

The instanton part of W is given by the critical value of the following functional: 

W inst (a;q) = 



(6.1) Crit^ ~ / p(x)p(y)G(x - y) + [ \p{x)tp(x) + Li 2 (qQ{ 

z Jexe Je 1 v 



'exe 
Equivalently: 

1 



<p(x)log[l-qQ(x)e-^) + Li 2 [qQ(x)e~^ 



-<p(y) 



(6.2) W inst (a;q)= f 
where ip{x) solves a TBA-like equation: 

(6.3) ip(x) = J G(x-y)log(l-qQ(y)e 

The solution is determined by the choice of the contour C, and the functions G(x), 
Q(x). Given this data, it is straightforward to solve (|6.3[) recursively: 

CO CO 

(6.4) <p{x) = ^ qVfa), P(x) = ^2 q k pk(x) 

k=l k=l 

6.2. The examples. We now give the expressions for the functions Q(x),G(x) and 
the contour C, for our examples. 

6.2.1. Periodic Toda. For the periodic Toda chain one has a simple formalism with 

r>( \ 1 ni \ d (x-e) 

P{x)P(x + e) ax (x + e) 

1 / A \ N N 

(6.5) WP crt = ^log (-)£<£ + £ ^ - a ") 

^ 6 ' n=l Z,n=l 

and the contour C that goes around the points a; + fee, k > 0. It can be deformed to 
go along the real line, when ai £ R+ §, and line > 0. 

In the perturbative limit the Bethe equation derived from the superpotential 
assumes the form: 

1 V.i Oj-Oj ^ 



("> (7) ' -n 



r 1 



£ ) 
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Figure 1. The contour C 



which has the form of the ordinary Bethe ansatz equations for a system of interacting 
particles (x\, . . . , xjv) with the factorizable S-matrix, which coincides with that of the 
open Toda chain. The two-body S-matrix is that of the Liouville quantum mechanics. 

6.2.2. The elliptic Calogero-Moser system. Let P(x) = n^i( x — a «)> an d 

P(x — m)P(x + m + e) 



(6.7) 



Q(x) 



G(x) 



d 



log 



P{x)P{x + e) 
(x + m + e) (x — in) (x — e) 



dx (x — m — e) (x + m) (x + e) 
Then: W(a, e; q) = WP crt + W inst , 

A N N 

(6.8) W pcrt = — r ^ a 2 n + {^M - a n ) - w £ {ai - a n - m - e)) 



n=l 



Ln=l 



and W mst is given by (|6.2p with the contour C in the complex plane that comes from 
infinity, goes around the points o/ + ke, I = 1, . . . , N, k = 0, 1, 2, . . ., and goes back to 
infinity. It separates these points and the points ai + Im + ke, I G Z, k = —1, —2, .... 
The perturbative Bethe equations have the form similar to (|6,6p : 



(6.9) 



n 



-m—ai+aj 

s 



-m+ai—aj 
e 



which corresponds to the system of N interacting particles on the circle of the size 
oc logq, with the factorizable S'-matrix of the hyperbolic Sutherland model. 
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6.2.3. Ruijsenaars- Schneider model. Let us denote: q = e _/3e , t = e~@ m , z = e~@ x , wi 
e P a i. The perturbative contribution is: 

N N 

W p Cr t = ™ J- a l + J2 (Piivn/O - ^(twi/wn)) 

n=l l,n=l 

where 

1 w n \ -« /log(to) + 27Tin\ 



n >) = ^E^T3^)~E^( — ^ — J 

^ n=l v ^ 7 neZ v ^ 

The P,Q,G functions are given by: 

(6-io) p(*) = fici - = p( ^y 

G(x)= d^^Cz- *«)(!- te)(l- 9 z) 



cte \ {% — t){z — g)(l — qtz) 
We assume |g|, \t\, \wi\ < 1, and the contour C is the unit circle \z\ = 1. 

6.2.4. XTie spectrum of observables. The spectrum of the type A quantum system is 
given by the solutions of the equations (|2.17p with W(a,e) given by (|6.1|) . Let us 
denote the extremum of the functional (|6.1|) by paOe)i ^a(^)- The spectrum of the 
type B quantum system is given by the solutions of the "dual" equations 

(6.11) cij = erii , rii E Z 

With these a^'s one can again define the functional (|6.ip and study its extremum 
Pb(x),ipb(x). 

The eigenvalues of the quantized Hamiltonians 

N 



H k = j &zdzti<S>{z) k = Y^p k l +. 



l=i 

are given respectively, for the type A and type B models: 

N 

(6.12) £ fc = V af + k / dx((x + ef- 1 - a;*" 1 ) ^(i) , 

z=i 76 V ' 

In particular: 

(6.13) £ 2 =£q-^-W(a;q) 

aq 
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6.2.5. On the relation between the type A and the type B models. The K = 2* gauge 
theory becomes the N = 4 super- Yang- Mills theory in the ultraviolet. The latter 
has the celebrated Montonen-Olive 6X2 (Z) symmetry. In particular, the 5-duality 
transformation maps the gauge theory with the gauge group G and the coupling r to 
the gauge theory with the gauge group L G and the coupling -r-, for some integer h = 
1,2,3. When the theory is perturbed by the mass term, the ^-duality symmetry still 
acts. We claim it maps the type A model of the elliptic Calogero-Moser system with the 
modular parameter r to the type B model with the modular parameter —-. The special 
coordinates a, map to raj. The modularity of the effective twisted superpotential is 
clearly supported by the expansion (|6.17p . 

6.2.6. More details and the origin of integral equation. Here we give more details in 
regard to the origin of claims from previous sub-sections. We do it for the case of the 
N = 2* theory (the relativistic case is studied analogously, along the lines of [171 14"5] ). 
One starts with the contour integral representation [TJ [U [7J [6] for the instanton partition 
function in the general O-background: 

(6.14) Z""'(a;q,m, El ,£ 2 ) = 

til n ^MiiQ(M e{m+£ f m ^f ^:, 

k=o Jt ^ l<I<J<k 1=1 1 z K ' 



where (pjj 



- PJ, 



x 2 (x 2 — e 2 )(x 2 — (m + £i) 2 )(x 2 — (m + £2) 2 ) 



(x 2 — e\)(x 2 — £2)(x 2 — m 2 )(x 2 — (m + e) 2 ) 

and Q(x) was introduced in (|6.7p . Next use the observation, reported earlier in 
118] . that (|6.14p is a partition function of a one-dimensional non-ideal gas of particles 
(pi, . . . ,<ftk subject to the external potential 

(m + ei)(m + E2)e" 



(6.15) U exb (x) = -log(Q(x)- 

and a pair-wise interaction potential 

(6.16) Mnt(x) = -log(D(x)) 

The free energy of this gaz can be analyzed using Mayer expansion |51[ [52| [53] . The 
subtlety with the £2 - limit is the clustering of the instanton particles, which leads 
to the multiple vertices of a given valency, labelled by an arbitrary positive integer k 
(the number of instantons in a given cluster) weighted by the partition function of a 
simpler one-dimensional gas 




n 



n.i 



l<KJ<k <t>IJ + £2 

which is equal to p^— , cf. [5]. By summing over /c's one arrives at the dilogarithm 
function in (|6.2p . Finally, the limit in (|3.7p singles out the tree diagrams, which lead 
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to (16. 2p . (|6.3p with the propagator function G(x) given by: 

s t- ■ D(x)-1 
G(x) = Limit £2 _>o 

6.3. The sum over partitions. The U(N) instanton partition function can also be 
written as a sum over iV-tuples of partitions A^ 1 ) , . . . , , which represent various 
configurations of instantons sitting on top of each other in R 4 . This representation can 
be effectively used to compute the first few terms in the q-expansion of W inst , or the 
first terms in the expansion, say, in the region |oj — aj\ 3> |e|: 

-minst/ \ Nm(m + e) , . 
W mst (a,e;m,q) = y — -ip(q) + 



m 2 (m + e) 2 d(p(q) 



(6.17) +-V : — q—T— + ■ 

where 
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7. Discussion 



To conclude we briefly review the relation of our results to some recent work. 

The periodic Toda system: In [19] the wavefunctions of the periodic Toda chain 
are constructed in the form of integrals involving the solution Q(x) of the Baxter 
equation. In order for the wavefunction to belong to the L? Hilbert space the zeroes 
6i,...,5n of Q must obey the quantization condition (in [19] the solution of Baxter 

equation is divided by the periodic function YliLi sinh y—^))- We claim that the 

quantization conditions of [19] are equivalent to our Bethe equations for the pure N = 2 
U(N) gauge theory, with the identification Si = a%. We have therefore provided the 
YY function for the quantization conditons in the periodic Toda chain. Note that our 
perturbative equation (|6.6p is derived for the periodic Toda system in [56] from the 
approximate analysis of the Baxter equation. The observation of [56J was very helpful 
in relating our two dimensional story [3J [5] to the four dimensional one. 

In [57] a quasimap partial compactification of the moduli space of holomorphic maps 
of a sphere into the affine flag variety LG/T is studied, and the corresponding J- 
function [28j is shown to obey tautologically a non-stationary version of the periodic 
Toda equation, the affine analogue of the result [29] for the G/T type A sigma model. 
The intersection homology methods of [57] are not applicable, it seems, to other quan- 
tum integrable systems. The questions addressed in the current paper do not appear 
in [57], while the interesting setup of [57] has its own place (14.10p in our general story, 
having to do with the four dimensional gauge theory in the presence of defects [50], 
which we discuss elsewhere [13] . 

The elliptic Calogero-Moser system: In [20] the quantum iV-particle elliptic 
Calogero-Moser system for the A^-i system, for the integer coupling parameter u was 
studied using the critical level limit of the free field representation of the conformal 
blocks of the WZW conformal field theory on a torus. The wavefunctions in [20] are 
written in terms of 



parameters ta i: on = 1, . . . ,i(y — 1), i = 1, . . . ,N — 1, obeying the elliptic Bethe-like 
equations, which are derived from the Yang- Yang function 



with the quasimomenta £i, . . . , £jv being the fixed parameters. One can supplement the 
Bethe equations of [20] by the quantization conditions £i = $+rii, £2 = "#+^2, ■ • • > Cat = 
1? + njv, with «j S Z. We conjecture 



Analogously, the critical value at fixed £ gives, conjecturally, the Yang- Yang function 
for the type B model. Note that the potential S(t; £; r) of [20] corresponds, in the table 
of dualities sketched in [4j [5] to ( a limit of) the four dimensional quiver gauge theory 
with the gauge group U (v — 1) x U (2 {y — 1)) x . . . x U ((N — l)(v — 1)), compactified 




(7.1) 



Crit (i ) V i(.-i) m S{t^;r) = I^ eff (ai,. . .,a N ; -ue,£,r) 
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on the elliptic curve E T . The equivalence 117. ip suggests an interesting duality between 
the 3\f = 2* theory in the special background with the coupling r and the quiver 
gauge theory. 

In [54] a partial resummation of the quantum mechanical perturbation theory around 
the trigonometric Sutherland model is proposed, leading to an equation on the eigen- 
value of the Hi operator of the elliptic Calogero-Moser system, for an arbitrary value 
of the coupling v. The relation of this approach to the YY- function formalism is not 
clear to us at the moment. 

Parallel developments: There are several recent developments which we feel 
are related to our story, including the obvious ones, such as the integrability in the 
AdS/CFT context, but also (43J EU EH Effl EE] . It is tempting to conclude from 
the convergence of these topics that the planar N = 4 super- Yang-Mills is equivalent to 
the union of the vacuum sectors of all gauge theories with eight supercharges (N = 2 
in four dimensions). 

We would like to return to all these questions in the future. 
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